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bl' Abstract 

We explore locally rotationally symmetric Bianchi I universe in 
Brans-Dicke gravity with self-interacting potential by using charged 
viscous cosmological string fluid. We use a relationship between the 
shear and expansion scalars and also take the power law for scalar 
field as well as self-interacting potential. It is found that the resulting 
universe model maintains its anisotropic nature at all times due to the 
proportionality relationship between expansion and shear scalars. The 
physical implications of this model are discussed by using different pa- 
rameters and their graphs. We conclude that this model corresponds 
to an accelerated expanding universe for particular values of the pa- 
rameters. 
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1 Introduction 

Current observational measurements obtained from many astronomical ex- 
periments (like Supernova (la), WMAP, SDSS, galactic cluster emission of X- 
rays, large scale structure and weak lensing etc.) strengthened the picture of 
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cosmic expansion at an accelerating rate PP-|1]. This significant phenomenon 
of cosmic expansion is prompted by dark energy (DE), a mysterious unusual 
kind of matter containing negative pressure. This is inconsistent with the 
strong energy condition and plays a dominant role in the composition of our 
universe [5]. The investigation of its obscure nature is one of the most fas- 
cinating issues in modern cosmology. Consequently, enormous DE proposals 
including Chaplygin gas, quintom, k-essence, phantom, quintessence, cosmo- 
logical constant etc. have been suggested [6l [7]. However, none of them 
provides an unambiguous solution to this problem and thus leaving it as a 
mystery for the researchers. 

General Relativity, in spite of its success in many ways, remained unsuc- 
cessful for describing the reality of DE and some other cosmological issues. 
This suggested the exploration of alternative theories of gravity by taking 
some modifications in Einstein Hilbert action [8l|9]. For this purpose, nu- 
merous modified gravity theories like Gauss-Bonnet theory, scalar tensor the- 
ories, f{R) and f{T) gravities and recently, f{R,T) gravity theory etc. have 
been constructed [101 E]. Among these modified theories, scalar-tensor the- 
ories are the most viable and interesting candidates of DE. In scalar tensor 
theories, the gravity effects are discussed by a tensor as well as a scalar field 
[12]- [13]. There are different scalar tensor theories available in literature, 
including the scalar tensor theories formulated by Brans and Dicke, Lyra, 
Nordtvedt and Wagoner, Saez and Ballester which are of particular interest 

Brans-Dicke (BD) gravitational theory is the most prominent and prevail- 
ing case of scalar-tensor theories which provides convenient solutions to many 
cosmological issues such as universe inflation and its late time behavior, co- 
incidence problem, cosmic acceleration etc [121 [SO]- Dynamical gravitational 
constant {G = ^), non-minimal interaction of scalar field with geometry, 
compatibility with the Dirac's large number and Mach's hypotheses as well 
as weak equivalence principle are some major facets of this theory [T2]-|14]. 
Various versions of BD theory are available in literature like generalized and 
chameleonic BD theory etc. for different cosmological implications. Brans- 
Dicke formalism attracted many researchers to explore exact cosmological 
universe models for cosmic expansion [21]- [23]. Recently, we have investi- 
gated exact solution of the BD field equations by using perfect, anisotropic 
and magnetized anisotropic fluids as matter contents [24] . 

According to grand unified theories, the phase transition (caused by the 
reduction in the temperature below some critical temperature in the initial 
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epochs of the universe evolution) results in the creation of some topologi- 
cal defects including domain walls and cosmic strings etc. These defects are 
responsible for density fluctuations and hence lead to a precise description 
of structure formation [25j. Since the cosmic strings interact with gravity, 
therefore it would be worthwhile to canvas the string's gravitational and as- 
trophysical consequences. The study of magnetic field effects in the matter 
distribution is of considerable interest as it provides an effective way to un- 
derstand the initial phases of cosmic evolution. Bulk viscous effects in the 
fluid lead to negative energy field and hence have a significant impact on the 
dynamics of the universe [261 [27]. Thus the study of magnetized cosmic stings 
under the influence of bulk viscosity leads to a better understanding of the 
dynamics of the universe. Various string cosmological universe models have 
been investigated in general relativity and scalar tensor theories |28j-|30]. 

This paper deals with the exact BD universe model with magnetized 
viscous cosmic strings. The paper is designed in the following layout. Next 
section provides BD formulation in the presence of self-interacting potential 
for Locally Rotationally Symmetric (LRS) Bianchi type I (BI) universe model 
and bulk viscous clouds of strings with electromagnetic effects. In section 3, 
we construct the universe model by solving the field equations. We discuss 
various physical parameters for the universe model. Lastly, we present a 
summary of the obtained results. 



2 LRS Bianchi I Model and BD Formulation 



In a simple BD theory, the BD coupling parameter remains as a constant, 
while its modified versions can be obtained by introducing variable BD pa- 
rameter, i.e., uj{(f)) and a self-interacting potential term. In Jordon frame, 
the action for self-interacting BD theory [23] is specified by 



Here is the matter contribution and L^j, is the Lagrangian density having 
BD scalar field as source and is given by 



where u is the BD coupling parameter (which is taken to be constant), R 
is the Ricci scalar and U{(j)) represents self-interacting potential. The self- 
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interacting BD equations obtained by varying the action with respect to 
scalar and tensor fields and are given by 

Equation (jlj) provides the evolution of scalar field. Here T and □ represent 
the trace of the energy-momentum tensor and the de'Alembertian operator, 
respectively. This theory leads to other modified theories when BD coupling 
constant takes some particular values |311 132]. In the limit, = 0o, a; — t- oo 
and T 7^ 0, the respective action and hence the field equations of GR could 
be recovered [33l 

In order to investigate the universe formation and the initial epochs of 
cosmic expansion, the study of Bianchi universe models is of great significance 
[35| |36] . The LRS BI universe model, as the simplest generalization of FRW 
universe model, is described by an anisotropic spacetime exhibiting spatial 
homogeneity [37] 

ds^ = dt" - A\t)dx^ - B^{t){dy^ + dz^). (5) 

Here the expansion in x direction is measured by the scale factor A, while in 
y and z directions it is measured by the scale factor B. The corresponding 
field equations ([3]) become 

2Ab B^ a ^5,0 a; 02 Uicj)) Too 

2^ + ^ = -2^t _ _ ^ , _ ^ 

B B^ 5 2 02 20 ' ^ ' 



B A Ab a bJ) uj(j)'^ 4> U{(t)) Tr 



B ' A^ AB ^A^ b\ 2 02 ^ 20 ' 
The BD scalar wave equation yields 

A B-^ T (20?7(0) -0m 

+ 3 - + 2- 0= ^ ^ — (9) 

^ M b'^ 2a; + 3 (2a; + 3)0 
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It is interesting to mention here that the energy conservation, T.^^ = 0, leads 
to a hnearly dependent equation as it is an outcome of covariant divergence 
of the BD equations ([3]) and (jl]). Thus we leave it and take the field equations 
©-© only. 

The average scale factor, the mean and directional Hubble parameters 

are 

a{t) = {ABy/^ H{t)=^-{j + 2^), (10) 

Hx = ^, Hy = H^ = —. (11) 

The anisotropy measure of expansion A, volume V, expansion scalar 0, 
deceleration parameter q as well as shear scalar cr can be written as 

3 

^ = V = a%t) = AB^ (12) 

i=l 

^ A B d A , I A B, , , 

e = < = ^ + ^B^ ' = S'77'-'- "^"TI'i-b'- 

Here A = corresponds to isotropic expansion of the universe model. 



3 Model for the Magnetized Viscous Cosmic 
String Fluid 

In this section, we discuss BI model with electromagnetic bulk viscous cloud 
of strings as background fluid distribution given by the energy-momentum 
tensor [27] 

t; = {p + Peff)u-^u, - p^fjs; - xx^ + e;:, (14) 

where is the particle's four velocity, A is the string tension density, E^^ 
denotes the electromagnetic part of energy-momentum tensor and is the 
spacelike unit vector that provides the string direction. Here we take = 
(1, 0, 0, 0) and = (0, 0, 0) which satisfy the relations 

u^u^ = 1 = -x'^Xf,, u'^Xf, = 0. 

The effective pressure Pe// is defined as the sum of isotropic and viscous 
pressures, i.e., Pe// = Pi + Pvis- We consider here the dust case, i.e., P/ = 0. 
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Moreover, Pyis = —^0, where ^ denotes the bulk viscosity coefficient. The 
electromagnetic part of is given by 

e; = nmux + 1/25^ - Kh% (15) 

where /i^ is the magnetic flux vector 

h^. = ^e^^^^F^^u-", (16) 

the terms JI, F"^ and efj^i^ap denote the magnetic permeability, the electro- 
magnetic fleld tensor and the Levi-Civita tensor, respectively. Moreover, p 
represents the proper density of strings being the sum of the particle density 
Pp (as the particles are attached to these strings) and string tension density 
A is specified hj p = pp + X. 

We assume that the magnetic field is generated in yz plane as its source 
is the electric current that fiows in x direction. Here the magnetic fiux vector 
has only one non-zero component hi. Moreover, the assumption of infinitely 
large conductivity along with finite current (in magnetohydrodynamics limit) 
leads to Fio = -F20 = -F30 = [38]. Thus, all the electromagnetic field tensor 
components vanish except F23. Using Maxwell's equations 

we find -F23 = n, a constant. The matter tensor (fT^ has the trace T = 
p+A— 3^G. The non-zero component of magnetic fiux vector will he hi = 
yielding the corresponding components of E'' 



-n 



2 



3 
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There are seven unknowns namely 0, U, A, B, p, A, ^ and only three inde- 
pendent field equations. For a closed set of equations, we take the following 
assumptions: 

• A = 5"; n ^ 1, this condition is constructed by setting the ratio ^ 
as constant, i.e., we assume the proportional relationship between the 
expansion and shear scalar |39]-| 



(p = 4>oB , a well-known power law relationship which indicates that 
the evolution of BD scalar field (p is dependent on the scale factor. 
Moreover, 6 > for expanding solutions 
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• U = f/o0^, a power law ansatz for potential in which k is any non-zero 
integer which can further be written as f/ = UqB^^ . 

• Also, we take = Mi, by setting the expansion scalar inversely pro- 
portional to bulk viscosity coefficient [27] which means that the rate of 
cosmic expansion decreases as the viscosity increases. Here Mi is any 
positive constant. 

The field equations for the fiuid ( IT^ turn out to be 

, ^Bd) ud)"^ UU) p n2 



E 52 B(f) 2 02 ' 20 2/1540' 

(18) 

5 2^2 50 a;02 £0 

\- n — - = — n + 1 

B 52 ^ 'b 6 2 6^ 6 26 6 



where we have used A = 5". 

Now there are three independent field equations and three unknowns 
namely p, B and A. Equation (fT9|) yields 

(2n(n + (5) + (wo + 2)5^)^2 2MiB^'-^^ u^b-^^+^^ 
25 H ; — -— = ——rz -T + 



{n + 6 + l)B 0o((5 + n + l) 7I0o(5 + n + l) 



+ 



MS + n + 1) 



whose solution is 
5 



2 2Mi5(2-5) n25-(^+2) 



0o((5 + n + 1)(2 - 5 + s) 7I0o((5 + n + + 2 - s) 
[/o^(^fc+2-^) M2 

^ 0o((5 + n + + 2 - 5 + s) ^ 5^' ^ ^ 
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where M2 is a constant of integration (taken to be positive) and s is a constant 
given by 

2n{n + 5) + {uj + 2)5'^ 



s 



(n + (5 + 1) 
Equation (pT]) leads to 

2MiS(2-5) n25-(^+2) 



dt 



'(f)o{6 + n + 1){2 - 6 + s) JlM^ + n + l){6 + 2 - s) 

^ 0o(5 + n + l){6k + 2-6 + s)^^^ • ^^^^ 

Introducing the notions B = r for time and X, F, Z for the space 
coordinates, the corresponding BI universe model becomes 

as -- 



>o(5 + n + 1)(2 - 5 + s) 7I0o(5 + n + + 2 - s) 

(/)o(5 + n + l)(5A; + 2-a + s) r'* ^ 
- T^dY^ + dZ^). (23) 

Equation (|T71) yields the density as 

p{B) = [(1 + 2n) + 5{n + 2) - '^]^,B^ B^'-'^ ^ 



2 2 27154 ■ 

Substituting the value from Eq. fl2T]) . it follows that 



p(r) = 0o(l + 2n + 5(n + 2)-^)-^ + [(l + 2n + 5(n + 2) 

io6^ 2Mi (l + 2n + ,5(n + 2)-^) 

~\5 + n + l)(2-5 + s)^ ^ (5 + n + l)(2 + 5-s) 

2^7Ir4 ^2 (5 + n + l)(5A; + 2 + s-5)^ ° " ^ ' 

The string tension density A can be expressed from Eqs.( fT8il and ( l20i) as 



(4 + 2w) ' 2 r2 0o(3 + 2u;) 

Mi(6 + 2w) ^ p 

X (l + 2uj + k)- , ^ , ^ + ^- — , 25 

^ ^ (3 + 2^)00^-^ 27Ir^+Vo 4 + 2w' ^ ^ 
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where 

f _ 1 Mi{2-6)t-^ (5 + 2)n2r-('5+4) 

T ~ {S + n + iy 0o(s + 2-6) ~ 2jI(f)o{s -6-2) 
Uo{6{k-l) + 2)T^(''-'\ 



+ 



20o(5(A; - 1) + 2 + s) 



p is given by Eq. fl24|) and ^ can be calculated from Eq.f l2T]) . 
The density of the particles is 

p{5 + 2uj) {3 + 2uj)<f)or\2r , ,^ u6^ 



A + 2uj A + 2uj T ' 2 



T 



2 



(3 + 2u;)0or^ ^W^-^^^o 
The directional Hubble, mean Hubble and deceleration parameters become 

H^: = nHy = nH^, = n\ ^ ( + + ) 



H 



bo{6 + n + l){2-6 + s) Jl(t)Q{6 + 2-s) 

(t)o{6 + n + l){k6-6 + s) Ti'+^)^ ' 
(n + 2) 2MiT-^ H2(5 + n + l)-i 



3 >o(^ + n + l)(2-5 + s) Jl(t)o{6 + 2 - s) 

(t)o{6 + n + l){k6-6 + s) Ti'+^)^ ' 

3(g + l) ^ 3r^ Ml HV-(^+^) 

+ — • 

200 ^ 

The expansion scalar is 

2Mir-^ H2r-(^+4) 



e = (n + 2)[ 



0o(5 + n + 1)(2 - 5 + s) 7I0o((5 + n + + 2 - s) 



M6 + n + l){k6 - 6 + s) T^^+^y ■ ^ ' 
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The coefficient of bulk viscosity takes tfie form 



Ml 2MiT-^ nV-(^+^) 

^ ~ (n + 2) ^(j)o{6 + n + 1)(2 - (5 + s) " JlM^ + n + l){S + 2-s) 

^ MS + n + l){k - 5 + s)^ ^(^^ ■ ^ ' 

The shear scalar is 

(1 - n) 2MiT-^ nV-('5+4) 



V3 >o(5 + ?^ + l)(2-5 + s) ~ 7I0o((5 + n+l)((5 + 2-s) 

(29) 



(7 = 



The volume for the model takes the form V — r""*"^ which is zero initially 

and becomes divergent when r — ?> 00 indicating that the expansion of the 
universe starts from zero to infinite volume. The anisotropic measure of 
expansion leads to A = ^(^^'^2)2 which is constant (as we have taken the shear 
scalar proportional to expansion scalar) and it becomes zero for n = 1. In 
our case, n 7^ 1, thus the universe model remains anisotropic throughout the 
cosmic time. The energy density remains positive for the allowed range of 
parameters as shown in Figure 1(a). However, it becomes infinite at initial 
epoch but decreases for final phases of the universe evolution. The energy 
density approaches to a positive value due to viscosity, given by 

■[l + 2n- — + (5(n + 2)] 



(5 + n + l)(s + 2-5) 



as T — )■ 00. If we take viscosity to be neghgible, then the density approaches 
to zero for later times thus representing an empty universe in future. 

The string tension density A increases with the passage of time and then 
remains a constant for future evolution of the universe as shown in Figure 
1(b). The particle density pp exhibits a similar behavior as does the energy 
density. It decreases from infinite value (at r = 0) to a certain constant due to 
viscosity effects and then remains constant for the future evolution as shown 
in Figure 2(a). The model results in dynamical deceleration parameter and 
consequently leads to negative as well as positive values for certain choices 
of parameters. For later time with k < 0, — 2 < u < —1.5, < 6 < 1 and 
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Figure 1: Plots (a) and (b) show the energy density p and the string tension 
density A versus time r. Here, red, green and blue lines correspond to u = 
—1.7, —1.8 and —1.9 respectively. Also, we have taken 5 = 2, Mi = 4, M2 = 
3, n = 3, fi = 0.5, k = —4, 11^ = 3, Uo = 1, and 0o = 1- 
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Figure 2: Plot (a) shows the particle density Pp versus time r and (b) indi- 
cates the deceleration parameter q versus BD parameter u. In graph (b), we 
have taken n = 2, the blue, green and red lines correspond to 6 = 0.2, 0.4 
and 0.6 respectively. 



n > 2, the deceleration parameter becomes negative, indicating accelerated 
expanding behavior of the universe model given by 

_ 3(g + l) _ ^ _ 3{2-6 + s) 
^~ n + 2 (n + 2)((5 + n+ 1)' 

However, for other choices of parameters, it exhibits decelerating behavior. 

The graph of deceleration parameter for particular choices of the parame- 
ters is given in Figure 2(b). It is clear that the model represents accelerated 
expanding universe as the deceleration parameter attains small negative val- 
ues as shown in this figure. The graphical illustration for expansion scalar 
is given in Figure 3(a), while the shear scalar, directional and mean Hubble 
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Figure 3: Plots (a) and (b) represent the expansion parameter 9 and viscosity 
^ respectively. Here, Mi = 4, Ms = 3, n = 3, /i = 0.5, k = -4, = 
3, Uq = 1, and 0o = 1- The red, green and blue colors correspond to 
0; = — 1.7, — 1.8 and —1.9 respectively. 

parameters exhibit a similar behavior. These parameters go to zero when 
r — )■ oo and diverge for initial time indicating the beginning of the universe 
model with a big bang explosion as shown in Figure 3(a). Figure 3(b) in- 
dicates the graphical behavior of viscosity parameter C, which is opposite to 
the expansion scalar (in accordance with our assumption). The viscosity pa- 
rameter is negligible at the initial epoch and increases with the passage of 
time as shown in Figure 3(b). Hence it prevents the universe to be empty 
for its future evolution. 

4 Summary 

It is argued that magnetic field has cosmological origin as there was highly 
ionized matter that was coupled to magnetic field which further leads to neu- 
tral matter as a consequence of universe expansion. It is interesting to discuss 
the magnetic field effects on the expansion history of the universe. Moreover, 
bulk viscosity effects and scalar field have important role for obtaining ac- 
celerated expanding universe model. Thus it would be worthwhile to discuss 
BI universe filled with magnetized bulk viscous strings for the discussion of 
early and late stages of the universe in scalar tensor theories. This paper 
investigates the cosmological model in self-interacting BD gravity with mag- 
netized bulk viscous cloud of strings by taking certain physical conditions. 
All the cosmological parameters depend upon the values of BD parameter 
u, parameter n as well as parameter 6 (that appear due to scalar field). 
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We have discussed the resulting model using different physical parameters 
through graphs. The results are summarized as follows. 

• According to Hubble 03], there would have been an infinitely hot and 
dense universe in its early phase. In fact, all the density of the universe 
was concentrated at a single point and hence initially there was zero 
volume of the universe. Penrose and Hawking [Il[|l5] have argued that 
expansion of the universe has been started from this dense and hot 
phase by an explosion and then universe is going to expand till today. 
After infinite time, universe would have infinite volume with negligible 
density. In our case, the proper energy density, string tension density 
as well as particle density remain positive for increasing BD parameter 
values and cosmic time. At initial epoch, the proper energy and particle 
densities diverge while they turn out to be finite for later time due 
to bulk viscous effects. Hence the presence of viscosity prevents the 
universe to be empty in its future evolution. Clearly, the physical 
behavior of our constructed model supports these arguments as energy 
densities turn out to be infinite (divergent) at initial epoch and hence is 
of considerable interest. However, the string tension density increases 
with the increase in time and BD parameter values. 

• Different parameters like if, Hx,Hy, O and a exhibit increasing behav- 
ior for decreasing r and become divergent at initial epoch indicating 
the big bang start of the model. However, these parameters approach 
to zero asymptotically, i.e., r — )■ oo. Thus the constructed model has 
an initial singularity. Moreover, scalar field evolves from zero (initially) 
to infinite value as r — )• oo. 

• As we have discussed in our previous work [23] , in case of BI universe 
filled with perfect fluid, there is a constant and positive deceleration 
parameter yielding decelerated expanding universe model. This prob- 
lem is then resolved by taking anisotropic matter contents. In present 
work, the deceleration parameter turns out to be a dynamical quantity 
rather than a constant due to the presence of self-interacting potential, 
viscosity and magnetic field effects which can provide negative values 
and hence yields accelerated expansion of universe. For example, if 
we take r — )■ oo and —2<uj< —3/2 as well as decreasing values of 
scalar field, the deceleration parameter gives negative values lying in 
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the range — 1 < g < which is in good agreement with the observed 
range for cosmic expansion as shown in Figure 2(b). 

• The anisotropic measure of expansion of the universe model turns out 
to be a constant showing that the universe model exhibits anisotropic 
behavior through the whole range of cosmic time. In [16], it has been 
pointed out that some large-angle anomalies are seen in CMB radia- 
tions, violating the statistical isotropy of the observable universe. For 
a better description of these anomalies, plane symmetric and homo- 
geneous but anisotropic universe models play a very significant role. 
Moreover, it is found [1T]-|19] that removing a Bianchi component 
in WMAP data can explain various large-angle anomalies yielding an 
isotropic universe. Thus the universe may have accomplished a slight 
anisotropic geometry in cosmological models irrespective of inflation. 
If the anisotropic nature of the model is maintained through the whole 
range of time, then it may explain or discuss these large-angle anoma- 
lies in CMB radiations. 

• The viscosity parameter increases with the passage of time which cor- 
responds to a non-empty universe in future. 

• In all physical parameters, the component of magnetic field has a neg- 
ative contribution, i.e., these physical parameters reduce due to the 
presence of magnetic field. 

In [no] , a class of Bianchi II, VII and IX filled with cosmic string fluid within 
Saez and Ballester gravity is discussed. It is found that the constructed 
model has no initial singularity which is inconsistent with the big bang model. 
Likewise, in a new work [51], the same class is discussed with cosmological 
strings in BD gravity and it is observed that the volume of the universe 
model is contracting rather than expanding which is physically unacceptable. 
However, in the present work, no such ambiguity exists and the results are in 
well agreement with the observations. It is worthwhile to mention here that 
our results are consistent with those already available in the context of GR 
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